History-dependent synchronization in the Burridge-Knopoff model 
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A three-blocks Burridge-Knopoff model is investigated. The dimension- 
less velocity- dependent friction force F(v) oc (1 + av) is linearized around 
a = 0. In this way, the model is transformed into a six-dimensional mapping 
x (^n) * x (in+i)j where t n are time moments when a block starts to move or 
stops. Between these moments, the equations of motion are integrable. For 
a < 0.1, the motion is quasiperiodic or periodic, depending on the initial condi- 
tions. For the periodic solution, we observe a synchronization of the motion of 
the lateral blocks. For a > 0.1, the motion becomes chaotic. These results are 
true for the linearized mapping, linearized numerical and non-linearized numer- 
ical solutions. 
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1 Introduction 

We consider the Burridge-Knopoff model (BKM) of three blocks p^. This is 
a simple chain of blocks connected with springs, with additional springs which 
drive each block with a constant velocity V. The blocks move on a surface 
with friction. The model belongs to a rich branch of many-body dissipative 
systems: the train model the Feder and Feder model |3] and the Olami- 
Feder-Christensen model 0]. There are also some connection from this branch 
to the Frenkel-Kontorova model [5] . At the root of the whole Tree of many-body 
nonlinear systems there is the fameous Fermi-Pasta- Ulam model [H]. There is 
also a conservative branch, with the Toda chain [7j, ding-dong [S] and ding-a-ling 
|5] models and presumably some others as leaves. 

Most characteristical effect displayed by the dissipative branch of the Tree 
is the stick-slip dynamics: the blocks move or stop at various time moments, 
periodically for one block and perhaps chaotically for two or more blocks. Recent 
works of Persson [13II14| allow to treat BKM as a generic model of two surfaces, 
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Figure 1: The Burridge-Knopoff system considered in this work. Three equal 
masses m are drawn along a rough platform with constant velocity V. The 
masses are connected with springs of strength k. Similar springs are used to 
attach them to the driving mechanism. 



one sliding on another with friction. On the other hand, the stick-slip dynamics 
has been observed experimentally in nanoscopic systems |15| . In particular, this 
kind of motion turned out to be relevant for the behaviour of tips of the scanning 
probe microscopes. 

The origin of the stick-slip dynamics has been often discussed in terms of 
the state and the rate dependent model 1171 IT%| . There, the character of 
the velocity dependence of the friction force, F(v), is considered to be crucial. 
On the other hand, as it was demonstrated in [2], the stability of the uniform 
motion depends only on the sign of the drivative F'(v). 

Usually, the discussion of all the above listed models concentrates on two 
topics: their application to theory of earthquakes (which are at the origin of 
the problem JHj) and the effects of self-organized criticality Both topics 
overlap at the Richter-Gutenberg law, which is an example of a scaling relation 
which has been observed out of laboratory ^21- Other effects under discussion 
within the above models are chaos and synchronization pQ. All these topics 
are often considered for large systems, despite the fact that even small systems 
of this kind are not entirely known. The aim of this work is to show some 
unexpected properties of the Burridge-Knopoff system of three blocks. 

In its original version pp the BKM is nonlinear because of two reasons. 
First is that the friction force F acting on a block depends on its velocity v in a 
nonlinear way, i.e. F — F(0)/ (l+cv). Second is that the equations of motion are 
nondifferentiable at time moments when a block starts to move or stops. In our 
previous paper jlHj we have demonstrated that for two blocks and small values 
of the parameter c, the first cause of nonlinearity is not relevant. Namely, we 
substituted the nonlinear function (1 + cv)^ 1 by its linear approximation 1 — cv. 
This step is justified for small values of c, i. e. for strong spring constant k, 
large mass m and low friction force which only weakly depends on the velocity 
v. On the other hand, this step allows to integrate the equations of motion 
analytically. It is useful to introduce here a new variable a = cF(0)/y/km. It 
has been shown in Ref. \HJ\ that for a < 0.3, the results of the original and 
linearized versions are practically the same. 
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Figure 2: Velocity dependent function of friction force F(v). Solid line for 
nonlinear form, dashed line for linearized form. Dotted line is a guide for eye. 

The goal of the present work is twofold. First observation is that the above 
mentioned linearization does not change the results of BKM for three blocks, 
when the parameter a is less than about 0.1. Second observation is that the 
character of the motion can depend on the history of the system. Near a = 0.1 
we observe a transition from periodic or quasi-periodic behaviour for small a 
to chaos for larger a. The character of the motion for a < 0.1 does depend on 
the initial conditions. Changing only the initial position of one of the blocks, 
we switch from a quasi-periodic and non-synchronized motion to a periodic 
evolution with the synchronization of two lateral blocks. This switching from 
one kind of motion to another is the same for the linearized motion and for the 
fully nonlinear formulation. In the latter case, the only solution is numerical. 
Here we use the Runge-Kutta method of fourth order. 

The paper is organized as follows. In the next section we describe the lin- 
earization procedure. The outcoming equations are piecewisely linear and they 
are solved analytically. The solution is valid except the time moments, when the 
number of moving blocks varies. These time moments are found numerically In 
Section III we show the results obtained within the linearized model and within 
its initial formulation - the latter is worked out numerically. The last section is 
devoted to conclusions. 



2 The linearization 



The starting point for the linearization procedure is the set of three equations 
of motion for three blocks 

mil = k(X 2 - 2Xi + Vt) - F(cXi) 

mX 2 = k{Xi-4X 2 + X 3 + 2Vt)-F(cX 2 ) (1) 
mX 3 = k(X 2 - 2X 3 + Vt) - F(cX 3 ) 

where F(cX) = F(0)/(l+cX), Xj is the position of i-th block, Vt is the position 
of the driving mechanism. The driving spring constants are chosen to be fc, 2k 
and k for the first, second and third block respectively, and the spring constants 
of the springs between the blocks are all equal to k. Equation for i-th block 
is valid as long as it moves. When the block stops, the condition of moving it 
again is that the total spring force acting on the block must be greater than 
F(0). 

If according to we change spatial and time units as follows: r = cot, Uj = 
kXj/F(0), uj 2 — k/m we can write the equations of motion in dimensionlcss 
form: 

U x = U 2 - 2Ui + vt - $(af>i) 

U 2 = Ui - 4U 2 + U 3 + 2vt - $(aU 2 ) (2) 
U 3 = U 2 ~ 2U 3 + vt- <P(aU 3 ) 

where v = V \Jkmj F(0) and a = cF(0)/^/km. The friction force has now the 
following shape: $(af/) = 1/(1 + all). 

For the linearized equations, the only modification is that the friction force 
$(a[/) = 1 — all. This approximation is justified as long as the constant a is 
small. Additionally, we have to check whether the friction force changes sign - 
this would be a harmful artifact. The profit is that the equations of motion are 
analytically integrable: the functions have this shape: 

U n = A exp(ar) + B exp(/3r) + C exp(7"r) + Dt + E (3) 

We need to solve a new set of equations whenever one of the block starts to 
move or stops. This time moments are roots of function (0 and they are to be 
found numerically. 

3 Results 

In Figs. 3a and 3b we show the time averages of the distance J(r) in the phase 
space between the first and third block. This distance is defined as 



J(r) = vW)-^(T)) 2 + (£/ 3 (T)-tMT)) 2 



(4) 
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Figure 3: The Euclidean distance J in phase space beetween blocks 1 and 3 as 
function of a. The solid line is for U\(t = 0) = +0.989 in the initial state, and 
the dotted line - for U\(t = 0) = —0.989 in the initial state, (a) for linearized 
equation of motion (analytical solution) (b) for nonlinearized equation of motion 
(numerical solution). 




Figure 4: The bifurcation diagram of W in the Poincare section in which W = 
(a) for Ux{t = 0) = -0.989 in the initial state (b) for U^t = 0) = +0.989 in 
the initial state. 



and its time average is performed over a time At = 7 x 10 3 after a transient of 
3 x 10 3 . In the numerical calculations, each time unit is 10 6 time steps. In both 
analytical and numerical calculations, one stick-slip cycle is about two or three 
time units. A comparison of the Figs. 3a and lb shows that there is a close 
correspondence between the results for the initial nonlinear formulation and the 
linearized one. In both cases, the synchronization appears for the parameter a 
between 0.005 and 0.08. In both cases, this synchronization does depend on the 
initial conditions in the same way. When the initial position of the first block 
J7i(t=0)=-0.989, i.e. it is close to the position for the stationary state, the 
motion is synchronized and < J(r) >= 0. (Here, < ... > is the time average.) 
This is so despite the fact that the stationary motion is unstable. However, for 
U\{t = 0) = +0.989, the motion is not synchronized and < J(r) > is positive. 

We note that the unstable character of the stationary motion can be proved 
in the same way as it was done for two blocks by Vieira 0. As long as the 
friction force decreases with velocity, some eigenvalues of the Jacobian are al- 
ways positive. The comparison of the Figs. 3a and 3b allows to believe that 
the linearization procedure does not alter the results for a < 0.1. Other re- 
sults presented in Figs. 4-7 are obtained for the linearized equations of motion 
and for the parameter a = 0.08. We compare the results for the two above 
values of the initial position of the first block. The initial positions of the 
other blocks are J7 2 (r=0)=-0.993 and J7 3 (r=0)=-0.991. The initial velocities 
are f7 x (T=0)=0.099, t/ 2 (r=0)=0.1012 and L> 3 (r=0)=0.1023. The initial values 
of the velocities are very close to the driving velocity v, which is set as 0.1. 

To construct the bifurcation diagrams, we use the sum of positions of the 
blocks with respect to the equilibrium positions: 

W = Ui - ut + u 2 - t/ 2 e + u 3 - Ui (5) 

We denote U{ = ur-5(l-ou)/6, Z7| = ur-2(l-ow)/3 and E7| = ur-5(l-ai>)/6 
the unstable point of equilibrium, which is found by taking U n = and U n = v 
in Eq. (J5J. The stability of this solution can be calculated as in In Figs. 
4a and 4b we show the bifurcation diagrams for the synchronized and non- 
synchronized motion. Subsequent points mean the Poincare surface section of 
W at W — 0. In the former case (Fig.4a, E/i(t=0)=-0.989) and below a = 0.1 we 
see a continuous line which means one fixed point of the mapping. This means 
that the whole motion is periodic. In the latter case (Fig. 4b, {7i(t=0)=+0.989), 
we see that below a = 0.1 the above line disappears. Instead, we observe a cloud 
of points for lower values of the velocity. This means the quasiperiodic or the 
chaotic motion. To check it, we have calculated the maximal Lyapunov exponent 
and we have found it to be close to zero. 

In Figs. 5a and 5b we show the time dependence of the distance J(t) between 
the lateral blocks in the two above cases. In the synchronized case (Fig. 5a) J(r) 
tends to zero as soon as for r = 4000, whereas for the nonsynchronized case J(r) 
displays stationary fluctuations around 0.14. Here again, the whole difference 
is the initial position of the first block. Finally, in Figs. 6 and 7 we show the 
character of the synchronized and non-synchronized motion. In the former case 
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Figure 5: The instantaneous distance in the phase space beetween the lateral 
blocks (a) for initial position of the 1st block U\(t = 0) = —0.989, (b) for initial 
position of the 1st block Ui(t = 0) = +0.989. 

(Fig. 6), the positions and the velocities of the block are periodic functions 
of time. In the latter case (Fig. 7), the motion is aperiodic. Moreover, the 
velocities of first and third block seem to be in some kind of antiphase: when 
one increases, another decreases. Both oscillate with the frequency more or less 
two times smaller, than the middle block. 

As we have seen in Figs. 6a and 7a, the stick-slip effect is present in both 
cases of periodic and aperiodic motion. The plots of velocities vs. time touch 
the value zero for each of the three blocks. It seems that some sticks are un- 
avoidable. This conclusion is supported by the fact, that the analytical solution 
for the linearized mapping contains the term exp[— F'(v)t], where F(v) is the 
velocity dependence of the friction force. Therefore, if only the derivative F'(v) 
is negative, the amplitude of the velocity oscillations increases with time and v 
must reach the zero value. 
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Figure 6: (a) Velocities of blocks as a function of time r. (b) Positions of blocks 
as a function of time t. The initial position of the first block U\(t = 0) = —0.989. 
Solid line for blocks 1 and 3, dashed line for block 2. 



4 Conclusions 

The linearization procedure, described in Section II, does not alter the system 
behaviour, if the value of the coefficient a is smaller than 0.1. This means, that 
the only relevant nonlinearity in this range of a is due to the stick-slip dynam- 
ics. Nonlinear character of the velocity dependence of the friction force is not 
essential for the kind of motion. This result is in accordance with our previous 
conclusions on the train model of two blocks [HJ\. From the computational point 
of view, such a conclusion is optimistic, because the time of calculations is much 
shorter for a map than for a set of differential equations. 

The history-dependent character of the motion below a = 0.1 can be inter- 
preted as an example of the frequency-locking, which precedes a transition to 
chaos [20] . When the initial values of the phases of the lateral blocks are close 
to each other, this frequency-locking leads to the synchronization. If they are 
oppose, we observe a kind of anti-synchronization, which persists in time. Sim- 
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Figure 7: (a) Velocities of blocks as a function of time t. (b) Positions of blocks 
as a function of time t. The initial position of the first block U\(t = 0) = +0.989. 
Solid line for block 1, dashed line for block 2 and dotted line for block 3. 

ilar effect are observed also in the Hamiltonian ding-a-ling system which can be 
described by coupled maps |2T). 

A new question appears, if the frequency-locking could be observed for a 
system with larger number of blocks, which imitates a real surface. This kind of 
behaviour is analogous to an excitation in linear systems. What are the conse- 
quences for the friction of macroscopic surfaces? At the beginning of this text, 
we have listed various many-body nonlinear systems. Could these excitations 
appear in some of this systems for larger number of subsystems? What are the 
necessary conditions for this kind of motion? To answer these questions, further 
research is needed, but the effort seems to be profitable. 
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